General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 




(NASA -CP-1532 20) THE DYNAMICS AND OPTIMAL N77-2 

CONTROL OF SPINNING SFACECPAFT WITH MOVABLE 
7SLE SCOPING APPENDAGES. PART C: EFFECT OF 

FLEXIBILITY DOPING BOOM DEPLOYMENT Final Oncla 

eport (Howard Univ.) 50 p HC A03/MF API G3/15 29086 


FINAL REPORT 


NASA GRANT 


Suppl. 1 




HOWARD UNIVERSITY 
SCHOOL OF ENGINEERING 
DEPARTMENT OF MECHANICAL ENGINEERING 
WASHINGTON, D. C. 20059 


FINAL REPORT 

NASA GRANT: NSG-1181 - Suppl. 1 

THE DYNAMICS AND OPTIMAL CONTROL OF SPINNING SPACECRAFT 
WITH MOVABLE TELESCOPING APPENDAGES 

Part C 

EFFECT OF FLEXIBILITY DURING BOOM DEPLOYMENT 

by 


Peter M. Bainum 

Professor of Aerospace Engineering 
Principal Investigator 

and 


Paul K. James 

Graduate Research Assistant 


May 1977 


ABSTRACT 


The dynamics of a spinning symmetrical spacecraft system 
during the deployment (or retraction) of flexible boom-type 
appendages is the subject of this investigation. The effect of 
flexibility during bona deployment is treated by modelling the 
deployable members as compound spherical pendula of varying 
length (according to a control law). The orientation of the 
flexible booms with respect to "die hub, is described by a sequence 
of two Euler angles. The boom members contain a flexural stiffness 
which can be related to an assumed effective restoring linear 
spring constant, and structural damping which effects the entire 
system. It is seen that the linearized equations of motion for 
this system, when Ihe boon length is constant, involve periodic 
coefficients with the frequency of the hub spin. A bounded trans- 
formation is found which converts this system into a kinematica lly 
equivalent one involving only constant coefficients. According to 
die Lyapunov reducibility theorem the stability of the reduced 
system can be analyzed using the standard techniques - i.e. an 
application of the Kelvin-Tait-Chetaev theorem. With the presence 
of structural boon damping only (and no damping on the hub) , 
criteria for stability: in the Lyapunov sense nay be developed and 
analyzed before or immediately after each deployment maneuver. The 
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NOMENCLATURE 


A 


b i 


C 


D 

6 



K 

k 


JL 

M, Mg, M3 
M 


dm(t) 


Linearized system state matrix 

Flexural coordinate in the transverse 
x,y plane 

Flexural coordinate out of the x,y plane 

Equivalent (structural) linear viscous rate 
damping of the booms 

Offset of booms from the principal x axes 

Synmetric positive senudefinite damping matrix 

Stew-synmetric matrix involving gyroscopic 
terms 

Composite principal moments of inertia 

<1, * i 2 > 

Hub principal moment of inertia 
(i = 1,2,3) 

Synmetric stiffness matrix 

Equivalent linear restoring spring constant 
of the booms 

Instantaneous boom length 

System composite parameters .in the equations 
of motion 

Synmetric mass inertia matrix 
Hub mass 
Tip mass 

Arbitrary differential mass within the 
system 


iv 


Stored mass 

Generalized force associated with the 
flexural coordinate in the transverse 
x,y plane (i = 1 4) 

Generalized force associated with the 
flexural coordinate out of the 
x,y plane (i * 1 + 4) 

Position of an arbitrary volume element 
along the boom (wire) 

Position vector of arbitrary volume element 
with respect to the hub center of mass 

Inertial position vector of the hub center 
of mass 

Inertial position vector of an arbitrary 
volume element within the system 

Offset of Idle booms from the system 
center of mass 

Bounded nonsingular transformation matrix 
used to transform nonautoncaxous system 
of equations to an autonomous set 

System kinetic energy 

Time 


Differential volume within the system 

Indicates the order of magnitude of ; the 
angular velocity vector components 

Naninal . value of hub spin rate 

Inerti tl coordinate system 

Bib coordinate system 

State vector 


5 

* 1 » e 2 
9 1 ,9 2 ,9 3 

P 

♦ 

Cm 

Superscripts 

rot 

Superscripts 

f 

H 

I 

♦ 

i 


Small perturbation associated with 
rotation about the spin axes 

Transformed form of the coordinates 
9 1 » 9 2 

Euler angles which define the orienta- 
tion of the hub with respect to the 
inertial system 

Boom (wire) mass per unit length 

Angle swept out in the transverse plane 
by the hub spin vector 

Inertial angular velocity vector of 
the hub 

Rayleigh dissipation function 

Indicates time differentiation 

s 

Associated with the final length of 
the boon 

Hub reference frame 

Inertial reference frame 

Refers to boon number (i = 1 4) 


j Index number of Lagrangian equation 

associated with a flexural coordinate 
Cj = 1, 2,3,4) 


I 

I 

I 

I 



CHAPTER I 
Introduction 

The effect of flexibility during boom deployment on the 
dynamics of a spinning spacecraft system is treated by modelling 
the deployable members as compound spherical pendula of varying 
length. The Euler angles, a^ and b^, describe flexural coordinates 
within the transverse plane and noraal to this plane (Fig. 1) . 

This model utilizing four deployable transverse booms (Fig. 1 * 
showing one boom) nominally lying in a plane orthogonal to the 
spin axes will be considered. The boons are assumed to emanate 
from points a distance, r, from the center of mass of the hub and, 
in their undeflected state, are not necessarily aligned with the 
hub principal axes (c^ ^ 0, 1*1 + 4). 

The development by Longman and Fedor^ for treating the 
dynamics of a system with constant length flexible wires will be 
modified to include the effect of varying the length. A quadratic 
approximation to the kinetic energy is developed which involves 
a total of 14 generalized coordinates. The degree of fle>dbility 
_is_ treated by introducing a restoring potential energy function 
proportional to the square of the flexural coordinates, where the 
proportionality (spring) constants can be varied according to the 
in-plane and out-of -plane resistance to bending. Boom bending in 
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the transverse plane and also out-of-plane can be simulated by 
varying the nagnitude of restoring linear spring constants , which 
can be related to the stiffness (El) of the boom material. 

A related recent paper by Janssens, considered the flexi- 
bility of fixed length appendages attached to a hub which was 
assumed to spin at a uniform rate with no transverse components of 
angular velocity. The motion of the appendages was simulated as 
spherical pendula in a centrifugal force field having two degrees 
of freedom, thus neglecting the perturbing effects of the wires on 
the motion of the central body. The author has examined this 
simpler problem to obtain the oscillatory behavior of the pendula, 
which describe the nature of non-linearities in the equations of 
motion, and gives a description of the asymmetry between forward 

2 

and backward swings in the presense of out-of-plane oscillations. 

Two other related problems have also been examined. The 
first considers stability boundaries on the extension of a pair 
of axial antennas whose undeformed state lies along the nominal 

3 

spin axis. It was assumed that the rate of extension was 
sufficiently small so that Coriolis effects due to the rate of 
change of length could be neglected. The author determines how 
far the antennas can be extended before a stability boundary is 

approached but does not simulate the actual dynamics during 

3 4 

deployment. In the second problem, the author's concern is to 


determine the maximum nutation angles expected to remain after the 


WWI'W ii " !r '>^ •: ■ :^'!!!Tl,l! : i^i^!' ; m! ! HlS!l'l!!'‘'"'^!: ! r«! 1 iW; , NW!1i^^^l^ ,;|! » ? ^ffll 11 - )! 


I 

1 

I 

I 


I 


f: 

li 

l 

l 

n 


r 

«; 


3 


deployment, the maximum bending moments, and deflections of long 
flexible booms. This analysis neglects boom extension by consider- 
ing boom leii£i~n (.%) to be constant in the formulation . The spin 

u 

rate (u» 3 ) is also assumed constant. 

The objective of the present investigation is to consider 
both the stability as well as the development of the deployment 
dynamics of spinning spacecraft with movable flexible appendages . 
The final objective will be to reach conclusions about the effect 
of first order appendage flexibility during extension or retraction 
maneuvers. 
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CHAPTER II 
Equations of Motion 

The first order equations of motion will be developed using 
the Lagrangian formulation, treating flexible boom (wire) members 
as spherical pendula of varying length. It is assumed that the 
amplitudes of the flexural coordinates as well as the amplitudes 
of the variational coordinates that describe the motion of the 
rigid part are small. It is sufficient to use a quadratic 
approximation tc the Lagrangian to obtain, the first order equations 
of notion. 1 We will first examine the development of the quadratic 
approximation to the kinetic energy. 

Fig. 1 illustrates an inertially fixed coordinate system: 

X, Y, Z, and a coordinate system: x, y, z, fixed to the symmetric 
hub of the spacecraft and centered about its center of mass. Also 
from Fig. 1, we let R^y be the vector from the cen • of the 
inertially fixed coordinates to an arbitrary volume element , 
dV, and R^, the vector from the center of the inertial coordinate 
system to the origin of the hub fixed axes. R^y, is the vector 
from the origin of the bib fixed axes to the arbitrary volume 
element dV. The quantity, dSgy/dt|j , is defined as the velocity 
vector of the volume element, dV, relative to inertial space. 


4 
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A. Kinetic Energy 

The kinetic energy of each element of volume can be expressed 
as, 

dT = dm|d^y/dt| 2 /2 (1) 

The total kinetic energy of tho. system is obtained by integrating 
Eq. (1) over the total volume in the system: 



where p is the density (mass per unit volume) . The hub and wires 
(boons) are assumed to be continuous, whereas the' tip mass and the 
stored mass are considered discrete particles. Thus, the kinetic 
energy for the entire system can be expressed as: 


T = (1/2) 


ilk 




dt I dt I 


dV 


hub 


* a '” /// "• ( 
wires ' 


/ I 


dt I dt I 


IdV 


, (1/7 , y / ^OV , • | \ 

tipLs^l « I 1 dt I 1 


/ HP 


+ (1/2) i dm- 
i=l 


dt |l * dt |I 
stored mass 


(3) 
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It should be noted that the last three terms in Eq. (3) include 
the effect of the varying length of the wires. 

B. Angular Velocity 

The orientation of the hub axes relative to inertial space 
can be expressed by using an Euler angle transformation , where 
0^, 0 2 , and © 3 are Idle Euler angles. The rotations follow the 
right hand sense, and are taken in the following order (1) a 
rotation through 0^ about the inertial X axes, followed by (2) a 
rotation 0 2 about the (intermediate) Y* axis, and then (3) a rota- 
tion 6 3 about Idle resulting Z" axis. After appropriate combination 
of the three transformation matrices , the following over-all trans- 
fornation describes the orientation of Idle hub system unit vectors 
relative to the inertial system unit vectors: 


V 


C0 2 C0 3 S0 3 C0 1 -S0 1 S0 2 S0 3 S0^S0 3 +C0^S0 2 C0 3 


X 



-C0 2 S0 3 cS^cSj+SO^sOjSOj -S0^C0 3 -C0^S0 2 S0 3 


Y 



S0 2 -sO^Oj c 9 1 c9 3 


_z_ 


where c, s indicate the cosine and sine functions, respectively. 
Fran. consideration of this specific Euler angle sequence. 


(4^0^+ 0 2 Y' + 0 3 Z” 


(5) 




^ ■■■ 'V ? v 'V. 
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Eq. (5) can be expanded by using Eq. (4) and the intermediate 
individual transformation natrices to yield the components of 


angular velocity) or, in the hub system as, 



r* ^ 

"L 


’ d i c9 2 c9 3 + 

0 2 S0 3 " 

_ L 



• 

• 

01 = 

“2 

~ 

“ 0 1 C0 2 S0 3 + 

®2 C0 3 


_“3_ 


t se 2 

+ ®3 . 


C. Development of First Order Expression for the Kinetic Energy 
In order to develop a first order expression for the 
kinetic energy, T, we can write: 

• 

0, = w(t) + <S0 O 
3 0 3 

e 3 = *(t) + <se 3 

where 50 3 /w Q and S0 3 « 1. In addition we assume 
l®ll> ! © 2 f a l a il J arjd l b iU <,: 1* For a rigid satellite in or 
close to its nominal state of spin (i.e. not during a deliberate 
spin-up maneuver) , 

w Q (t) = o* n = constant (about spin axis) 

(I 

<|>(t) = <* Q t . 

Since small angles have been assumed, we can replace, e.g. ■ 
Cos 0 2 by 1, sin © 2 by © 2 , etc. Then, with the aid of Eq. (7), 

Eq. (6) can be approximated by: 



( 9 ) 
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8^o|» + 02 s ^ 
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-0 1 «0 3 S<N-0 2 50 3 C1(! 
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• * 
-0jSl|> + OjOp 


• • 


-V 


- *®3 - 


®1 0 2 - 


W„ 


W, 


w„ 


where the three dots indicate terms of higher order than quadratic. 

The generalized coordinates for the center of mass of the 

• • 

hub are the cartesian components of in inertial space. X, Y, 

Z represent components in the coordinate system (I) of dRQg/dt|j, 

♦ • • 

an inertial derivative. The quantities Xg,yjj>and are defined 
as ex ponents of this vector projected into the hub (H) , coordinate 
system y^ z^) as seen from the inertial reference frame. 

The first term of Eq. (3), the kinetic energy of the hub, 
contains contributions due to the translational kinetic energy 
of the center of mss and the rotational kinetic energy about the 
center of mass. 


T hub = (1/2) ♦ *H + *H> + (1/2) + ”t> * (1/2) ^ 4 

( 10 ) 

It is assumed in the development of Eq. (10). that the hub is 
symmetric and homogenous (i.e. 1^ = 1^) . With the aid of 
Eq. (9), Eq, (10) my be expressed as, 


T hub = (1/2) + *H + V + (1/2) \ C( V Wi 

+ (Wq+W^) 2 ] 

+ (1/2) 1^ (Wq+W^) 2 . (11) 

As an example, Wq+W^+Wj may be expressed in terms of the Euler 

angles and rates as follows: 

. • • • • 

(W q +W 1 +W 2 ) 1 = 6, cos * + 0 2 sin * + (-e^e^in <|> + 0 2 50 3 cos v ) 

( 12 ! 


W, is higher order than quadratic; thus it can be emitted, along 
with all other, higher order terms. 

To calculate T wire and it is necessary to determine 

the inertial velocity, dR^y/ dt | ^ . From Fig. 1 , and with the 
application of Coriolis’ law, 


^dtli + + ffx 


( 13 ) 





(It should be noted that in the development of T^ tih this expression 
simplifies to: 

d^jy/dtJx = + S'. x ^ 

since dEgy/dt| H = 0 for a differential volune element inside the 
hub.) 

The quadratic approximation of R^y, written in matrix form, 
following Fig. 2, can be expressed as: 


[*Hv] H = ( ^ +r) 


- 



OQS C^ 


-a^ sin 

sin c^ 

♦q 

a^ cos 

0 


b i 

L J 

r~ 

’ “I 


-(l/2)q(b? + a?) 


COS C. 


I sin 


+ •••■ + (q*r3t 


cos c • 


sin c- 


-a^ sin c^ 


cos c^ 


b i 


+ 0/2) jja>? + a?) 


COS c^ 


sin c . 


+ ... 

where q = q(t) 


(15) 


In Eq. (15) , the ccnponents of R^y are given in the hub coordinate 
system, (H). 

For the case of dV along the wire at an instantaneous 


length, qCt) , 
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CdFjjy/dt = q 






™ **“ 

cos 


-a^ sin c^ 


-a^ sin c^ 

sin 

+ q 

4^ cos c^ 

+ q 

cos c^ 

0 


b. 


b, 



x _ 


L i J 


- q (b*b. + a*a*) 


COS Cj 
sin c . 


Cl/2) q (b? + a?) x 


cos c. 


sin c 


■i 

0 J 


(16) 


The quadratic approximation tc the kinetic energy per unit mass of 
the volume element can be obtained by calculating, 

(1/2) [dl^/dtl^ • [dS QV /dt | X ] H (17) 

through second order. 

If dV is along the wire, Eq. (13) may be expanded with the 
aid of Eqs. (15) and (16) to yield: 

<®Ov /dt >I = ~ qi i sin c i “ q (b i b i + a^) cos c^x 

+ Cy H + qa^ cos c^ - q(b^b^ + a. a.) sin c^]J 
+ Cq b i + Zjj]k + u x . 


( 18 ) 
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In order to include the effect of the wire on the composite 
system moments of inertia we consider, for example, the contribution 
of the wire to I y , 

V . - - * *» 2 (i9) 

■ ' ware 


where I , I and I represent the principal moments of inertia 
uu w ww r 

of the wire, (u is the longitudinal wire axis, w and w are transverse 

2 

axes). The term md represents the moment of inertia of the wire 
due to the offset of its center of mass from the hub principal 

A A A A 

"yjj" axis, o, B, y are the direction cosines: (u • y), (v *y), 

and (w • y) , respectively. Assuming the wire is thin and noting 
that for small displacements of a^, b^, y « 1, Eq. (19) may be 
approximated by: 



* I S 
w 


Z 


+ md 


( 20 ) 


2 2 

vhere 8 = cos a^ . The moment of inertia of the wire about the 

v axes (Fig. 2), can be expressed as: 

1/2 

1^ = 2 f _q 2 pdq = m£ 2 /12 (21) 

0 


where p is the mass per unit length. After substitution of Eq. (21) 
into Eq. (20), 


I * (nuT/12) cos a^ + m [r+(i/2)cos a.jT (22) 

wire 

which can be simplified to, 

I„ _ - (mi 2 /3) cos 2 a. + nr 2 + rart cos a, (23) 

y setre 1 i 


The contribution of the moment of inertia of the tip mass 
can be developed in a similar manner as: 

Xy - n^CCr^tcos a^) 2 + ( isin a^) 2 ] (24) 

y tip mass 

The combined effect of the appendages can be obtained by adding 
Eqs. (23) and (24) to yield: 

I - (m/3+m t )A 2 +2(m/2+m^)r*. + (m+m^r 2 (25) 

appendages 

for the case where « 1. From synroetry, the wire and tip mass 
in the (ix,-±;y) quadrant also contribute to I y . Eq. (25) now 
becomes, 

I y = [2(m/3+m t )£ 2 + 4<m/ 2+n^) Ar + 2(m t +m)r 2 ] = I y (26) 

Assuming the hub is symmetrical, (I„ = I„ ) then the composite 

H 1 H 2 

principal moments of inertia can be expressed as: 


*1 8 X 2 = \ + X y < 27) 

The composite moment of inertia about the third principal axis of 
the hub can be written as. 




+ I 


tip mass 


(28 



The second term in Eq. (28) carl be approximated as: 


I, + nd (29 

. wire -wire ' 

With the aid of Figure 2 , Eq. (29) can be expanded to yield, 

2 2 

I- = (at /3) + nr + mtr cos a* (30 

3 Wire • 1 

The third term of Eq. (28) can be written as: 

2 2 

I, ; . . = m^c + '-2a t ir cos a^ + t ^ (31 

tip aass 

After substitution of Eqs. (30) and (31) -into Eq. (29), 

Ij s Ijj^ + (a t ’*m/3)t^ + 2(a t +m/2)tr + (a t +m)r^ (32 

If we consider the four appendages, then Eq. (32) can be written as 
I 3 = Ij^ + C4(a t +m/3)t 2 + 8(a t 4 m/2)tr + 4(a t +m)r 2 ] (33 

When Eq. (17) is expanded in the calculation of T' w ^ reg , 
terms which are periodic such as sin^+c^) and cos(^+c^) appear. 
If we consider the four appendages, then the terms can be related 


using the following transformations : 
Let a ~ iji + c, thus, 



sin a = sin a 


cos a = cos a 

sin (a+ir/2) = CM a 

cos (a + n/2) s - sin a 

O' 

sin(a + id = - sin a 
cos (a + n) = - cos a 
sin(a + 3n/2) = - cos a 
• cos (a + 3ir/2) * sin a 

New, volume integration over each wire can be performed, 
with the result that the kinetic energy of the wires for constant 
length is calculated as: 

T . = (1/2)C2(0? + 0 2 )mA 2 /3 + 2(0? + 0?)nnZ 

wires 1 l i i 

+ 2(0 2 + OjJr^m] ( 

If dV were at the i— tip mass, then expansion of Eq. (17), 
multiplication by m t , and substitution of l for q (for constant 
length), would result in the kinetic energy for the i— tip mass. 
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*tip W s tl/2,( ^ * »H + *H> »t + W + (1/2)( *H ) V 2 

♦ (1/2)Wq + a?) - rA(b 2 + a?)] 

+ (l/2X«e 3 ) 2 n^Cr 2 + A 2 + rA) + w Q 9^ m t (r 2 +A 2 +2rA) 

+ w Q «e 3 m^r 2 * A 2 + 2rA) + (1/2) (a 2 Wq A 2 n^) 

• • * • • • 

- mj.ZCa^ + y H a^w Q ) sin + n^ACy^ a^ - Xy a^) cos 

» * • 

- n^Atb^Wg 9^(r+A) + b^ e^Cr+A)] cos C^+c^) 

m • • 

- n^ACb^Wg SjCr+A) + b^ 9^(r+A)] sin ('l'* 0 ^) 

+ n^Ce 2 + e 2 )(r 2 + 2rA + A 2 ) (36) 

If we consider variable length appendages, the fo rm ula t ion 
of Ref. 1 can still be used, except length is now a function of 
tine and the position of an arbitrary volume element along the 
wire, (q), is a function of time. The additional terms which 
result for dRgy/dt | deployment 0411 ** ex P ressed 33 follows : 

(deployment = C 9 =°s c. - 4 a. sin c. - 0/2) 4 <b? + a?) 

• • 

cos c^3 T + [q sin c. + q a. cos c^ 

- (1/2) q(b 2 + a 2 ) sin c. ] J ♦ q k 


( 37 ) 
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When Eq. (37) is substituted into Eq. (17), the resulting 
additional terms in the kinetic energy per unit mass are: 

'“daploywm a tl/2) * 4 “o « a i * « "t « «i 

* • • * • • * 

- q 0 ^ 2 ^ a i + a i a i + 3 w o a ita fp * + <1 a i 00 3 (q+r) 

• • • * • 

+ q Oi e^jCq+r) + 3^q cos c i - x^q sin 

- d/ 2 )Xjj q(bf + ap cos c^ + y H q sin cp cos c i 

- (l/ 2 )y^q(b? + a?) sin c^ + q + q q b^^ 

• * • • 

+ rqbp 2 "3 s ^ n ^ + c p + ^ s i^i sin(t|»+ cp 

• •, * • 

+ rq 6 ^ 50^ b^ cosCty + cp - r 0 2 q cos( cp 

(38) 

The limits of integration over the differential volume in the 
expression for T • _ can be formulated as follows: 


q(t) = qt, where q is assumed to be constant (for an 

* 

assumed uniform extension rate). Then, dq = q dt, and 


l l/q 

f(. . .)dq -*•/ (. . .)q dt, 0 < t < t. 

o J o J 


when 

when 


t = 0, q = 0; t = z/z, q = z(t) 


• • 


Z <_ Z^, t^ = zpz when t(0) = 0; q = Z 


where and represent final length and time, 
respectively. 

The following additional terms result for T . . 

^ wires 

. 4 # 4 

T - = (2/3)mZ 2 + (1/2) l w n a.m.4 2 + (1/2) Y a. 60-n. i.£ 2 

wxres depl. i=l 0 1 1 i=l 1 % 

4 4 

+ (1/2) J i a i d 1 0 2 rm i i + Cl/2) £ (39! 


If we consider T!Wgg , where the expansion of Eq. (17) 

is similar to that for then multiplication by m^ and substi- 

• • 

tuting i for q yields: 


tip mass 


depl„ 


• 2 4 ... ■ 4 . . 4 

= 2m J. + l tow n a.r + [ to a . 59 -r + £ 
z i=l ” u 1 i=l r 1 d i=l 


^ b i^ 

(4o; 


To consider the kinetic energy due to the stored mass of 
the rods inside the hub, the following assumptions are made: 

(1) at full extension there is no stored mass; 

(2) the mass is stored at a point inside the hub, 
a distance r from the hub center; 


(3) the extension rate is the same for all rods; 

(4) the wires are homogenous. 



The extended mass per unit volume along the wire can be defined 


m(t) = p£(t) (41) 

where 0 <_ t <_ t^, and p is the mass per unit length. The stored 

nass per unit volume can be defined as: 

dm. = (£ f - £(t))p (42) 

stored 

where £(t) resprsents the instantaneous length. Since the stored 
mas s is inside the hub, then Eq. (14) can be used to express 
dS^y/dtlj , where; 

- r(cos c^ l + sin c^ J) (43) 


If we substitute Eqs. (14) and (43) into Eq. (17), and apply the 
sane procedure, the terms which result for the combined effect 
of the four locations, ^ield: 


T 

stored mass 


= (1/2) .1 Pi^f^t)) [4(x^ + y 2 + z^) + 4r 2 Cw 2 

+ (0 1 0 2 > 2 + ( <$9 3 ) 2 } + 8r 2 (w Q 50 3 + Wg 0^02 + 

+ 8r 2 Q. ± 8 2 <S0 3 + 2r 2 (0 2 + 0 2 ) + 2r 2 (8^ + 0 2 )(60 3 ) 2 ] 


( 44 ) 


Eqs. (11), (35), (36), (39), (40) and (44) can now be 
combined to provide the quadratic approximation of the kinetic 
energy for the entire system; which can be written as: 



T = (l/2)M(Xg + y 2 + z^) + x 3 w Q ( 5 83 + S-^) + (Jl.) I3 w Q 4 

4 

+ (1/2) I-^e 2 + i\) + (1/2) I 3 («S0 3 ) 2 + ^(Wq + 683) J a i 

+ 0/2) M3 w 2 j i a 2 ♦ Mj % b. - (1/2) M 1 w 0 2 .^(a? ♦ bj 

+ (1/2) M- 7 (a. 2 ♦ b.) 

4 ^ i=l 1 1 . - 

+ Mj J C- ^(^ s:n c i + w Q a. L cos c i ) + y H (a i 


cos c . 


4 . 


- w Q a^ sin c^] + ^ J Cb^ sin (iji+c^ - Wq^ cos(*+c i )3 


• * 


- 9 2 ^i 003 ^ +c i^ + w q b i sin (t+c^)]) + (2/3)m l 


+ (1/2) £ w n a-rm- A 2 + (1/2) £ a. 66- rm.£ 
i=l 011 i=l 1 d 1 


:2 


4 


4 . 


+ (1/2) l a i 8 1 8 2 rm i l + (1/2) ^ 2 \ l 


4 


J x WV * a i d9 3 "V + Jj. Wl- "V 


4 . 


+ l + (1/2) l P i (Z f -£(t)) IXx 2 + + + 4r>4 Cw 


+ ( 0 1 0 2 ) 2 + (68 3 ) 2 ] + 8r 2 (w Q 683 + w Q 0^2 + 59 3 0 i 3 2 ) 


o • • • 0*9 *9 9*9 *9 9 

+ 8^0^2603 + 2r z (0f + 0 2 ) + 2r (6^ + Qp < 4 V ^ 


where the varying length terms (containing £) and the stored mass 
terms are additions to the constant length terms developed in Ref. 
1. In connection with Eq. (45) the following parameters 
have been defined: 

M = m^ + 4m + 

M^ = (n^ + m£/3 + n^r + mr/2)£ = Mg + rHj 
M 2 = (n^ + m /2 )£ 

Mg = (n^ + m/3 )£ 2 

where m = m(£) , and are defined similarly to those given in Ref. 1 
Terms associated with the stored mass in Eq. (45) can be 
included with the composite moments of inertia expressed in 
Eqs. (27) and (33), and also with M_, aiid Mg- These' 
quantities can now be redefined in modified form as: 

M = mg + 4p£(t) + 4m stored + 4m t = constant (46a! 

m stored = pU f - Z(t)) 

= (n^ £(t) + (l/3)p [£ (t )] 2 + nn t + (l/2)r£(t)) £(t) 

(46b! 
(46c: 



M 2 = (m t + (l/ 2 )p £ (t)) £(.t) 
M 3 * (ji^ + (1/3) p£(t)[£(t >] 2 


(46d 
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T l = ^ + C2(m t + (l/2)p£(t))[£(t)] 2 


+ 4 ( 11 ^ + (l/2)p£(t))r l (t) + 2 ( 11 ^ + pi^)r (46e) 


I 3 = Igg + C4C1H,. + (l/3)p£(t)C~(t)] 4 

+ 8 ( 11 ^ (l/2)p£(t))r£(t) + 4 ( 11 ^ + p£ f )r 2 ] (46f) 


For varying length, m has been defined in Eq. (41) , assuming the 
mass per unit length, (p), of each appendage is the same. 

The corresponding time rates of change of the (new) composite 


parameters and of the (hew) composite moments -6F ihertia ean be 


written as : 

M =0 

* * * 2 

= 2m^ ££ p ££ + m^.2. + rp££ 

• • • 

+ pll 

• . * n 

^ = 2m t U + p££ 

• - • *0 • • 

= 4m t ££ + 2p££ + 4m t £ + 4rp££ 

* - • *0 • • 

I 0 = 8m.££ + 4p££ + 8mjr£ + 8rp££ 

o t u 

• • 

where £ = £(t) and £ = £(t) 


(46g) 

(46h) 

(46i) 
(46 j) 
(46k) 
(461) 


D. Development of the Equations of Motion 

The fourteen generalized coordinates selected are; X, Y, Z, 

0 1* 0 2* 0 3* a l» a 2* a 3» a 4» b i» b 2* b 3 31x1 b 4* We note tbat X » Y 
and Z, are inertial ly fixed coordinates and do not appear in Eq. (45) , 
thus they are cyclic. Lagrange's equations of motion: 

d/dtOT/dq^) - 3T/ aq i = Q ± i = 1,2,... 14 (47) 

where far the wires (booms) is the generalized force derivable 
frcm a potential energy function which is assumed to be proportional 
to the square of the flexural coordinates, a^ and bp and involves 
a linear structural restoring spring constant, k; also Q is derivable 
frcm a Rayleigh dissapative function, '■f (i.e. "P oc a? and 
oc b^ , i=l - 4 bocms) , involving, assumed linear viscous 
structural rate damping, C. Eq. (47) nay be developed in the 
coordinates, q^, to yield: 

d/dt(3T/3X) = = p x ; d/dt(3T/3Y) - C 2 - p 2 ; 

d/dt(3T/3Z) = C 3 = p 3 (47a) 

where p , p^, p 3 are the constant generalized momenta associated 
with X, Y and Z, respectively. If we choose the inertial axes such 
that there is no momentum of the spacecraft relative to the inertial 
axes at time t = 0, then the constant momenta are zero, therefore 
Ci = C 2 = C 3 = Q. Eq. (47) can be developed for this system as a 


24 


set of fourteen equations of motion. Furthermore if we keep 
• • • 

Xg, yjj and as coordinates, then solving the linear differential 

equations with variable coefficients (when i is constant) , is made 

easier. X The equations of motion can be expressed as follows: 

4 . 4 „ 

“ ^t* + £ ^ a i sXn c i * w 0 038 c i^ “ ^2 I s ^ n c i 


+ w Q cos c i + w^ cos <^) =• 0- 

4 . 4 

y H + M2 l (a^ cos c^ - Wq a^ sin c^) + M2 £ (a^ cos c^ 

i" 1 1“ 1 


Wq a^ sin c,. - w Q sin c^) = 0 


4 


**2 X B i' + **2 . X b i + X - p i *- b i 11 + b i (£ * 
1-1 1—1 1—1 


ll) 3/2 


+ l m. (b • l + 1 b . ) = 0 
i=l ^i 1 1 


4 . 


h W 0 9 2 + X 3 (w 0 9 2 + w 0 9 2 } + I l 9 l + I l 9 l + ^. i ? 1 Cb i siz > + c i ) 

4 .. 2 

- Wq b^ cos(f + c i )] + M^ l [Cb^ + w Q b^) sin (# + c^) 


- Wq b. cos + c i >] + l rp ± (a 1 9 2 ££ 2 + a,. [9 2 

I'll 


# 2 

.9~ U 4 


+ 9 2 U 3 + 2£ ££)]}/2 + l rai . [a t 9 7 £ + a-(9 9 l +9 0 £)] 

i=l \ 1 2 12 2 


+ 4 J p i r 2 {z f (0 2 5e 3 + e 2 «e 3 ) - U 0 2 «e 3 + £(i 2 50 3 


+ 0 2 50 3 )}} = o 


4 


^1 0 2 + ®2 “ 003 ^ + c i^ + w o s * n ^ + c i^ 

* ~ 2 

-Mi I C(bi + Wq b^) cos (tjt + Ci ) + w Q b^ sin($ + c^)] 

4 

+ 4 .^ 1 P i r2{Jl f (0 l 60 3 + e i 69 3 ) “ ^ 9 i 50 3 + ^(9 1 59 3 +;9 i'50 3 )]> 

- I 3 Wq 0^^ - ( 1 / 2 ) l a^ © 1 rip 4=0 ( 


. 4 . 


[ 3 (w Q +60 3 ) + I 3 (w Q + «0 3 ) + M x J a. ♦ ^ J &i 


+ (1/2) l p.r(a. 44 2 + a. [4 d + 2 III]) + £ rm + (a .4 + a.i 

i=l i=l H x 1 


- 4 l p. ( 4 - -4)r 2 3,5, : 0 
i=l 1 r 12 


*l (w 0 a j + w o + 60 3 ) + ^L (w o + 6 ty + - w o a i } 


+ Mgaj + M 2 l CC-Xjj + y H w Q ) sin < 5 . + <y H + ^ w Q ) cos c^ 
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I 1 

I 

I 


r [ : l, ! 


+ Mg J (Xg sin + y H cos c^) - (1/2) J p .rW (w Q 
i*l i*l 

. . . . 4 

+ + 0^0 2 ) " 0 l®2 r 1 X “'t, + tej + s 0 (48g) 

.... . . 4 

MjZjj + M 2 Zg + Mjbj + Mgbj + £ CCe^ + 2w Q 0 2 ) sin(ij> + c^) 

- . . 4 

+ (- 0 2 + 2w Q 0^) cos (4» + c ± )] + l C9 1 sin ($ + c^) 


- © 2 cos (<|» + c^)] + Wq bj + kbj + Cbj = 0 


(48h) 


• • • A • 

where = 2n^ ll + p lir + rpil in Eq. (48g); 

• • 

Mj = CpJU/2) in Eq. C48h). 

In Eqs. (48g) and (48h), j = 1,2, 3, 4; i.e. these equations represent 
a set of eight Lagrangian equations corresponding to the coordinates 
a 1 - a^ and b^ - b 4 , respectively. The terms ka ^ , kbj, Ca^ and 
Cbj in these equations, are associated with the generalized force 
for the wires (boons; due to the effect of a linear (structural) 
restoring spring constant, k, and structural damping represented 
by C, (a linear viscous rate damping constant), which is assumed 


i 


to exist in the four boons. Both k and C are assumed to be the 
same for each bocnii and associated with a potential energy function 
and a Raleigh dissapative function, respectively. 


Since Eqs. (48a), (48b) and (48c) yield first integrals 
• • • 

where x^, and can be related to the constant generalized 

M M M 

momenta 'pp p 2 , and p 3 , respectively, then x^, y H and are zero 
and can be eliminated from these equations. Eq. (48) can be 


reduced to a set of eleven equations by solving Eqs. (48a), (48b) 

• • • 

and (48c) for x^, and z^ respectively and substituting than into 
Eqs. (48g) and (48h). In Eqs. (48), we will now assume that the 
tip mass, hl , and mass per unit length, p., are the same for the 

i 1 

four appendages. 
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CHAPTER HI 
Stability 


To examine stability before and immediately after deployment: 

• - i 

we assume l = l - 0. Gq. (48) for constant length , can be 

»• 

reduced to the following set of ten equations by solving for 66^ in 
Eq. (48f ) and substituting the results into Eq. (48g). 

I^l + ^3 m o® 2 + + w o sin + 

+ (b 2 + Wq b 2 ) cos (♦ ♦ <^) - (b 3 + Wq b 3 ) sin (i|> + c^ 


(b^ + b^) cos (4» + c^)] = 0 


(49< 


I l 0 2 * I 3 W 0 9 1 “ ^L^ b i + w o b 2^ 003 ^ + c i^ 


+ (b 2 + Wq b 2 ) sin(i|> + c^) + (b 3 + Wq b 3 ) cos (<i» + c^) 


(b^ + Wq b 4 ) sin C4» + c^ = 0 


(491 


4 . 


/I 3 )i. + (M 1 -M* 3 )Wq a^ + MgCC-Mg l (itj cos c . 


4 . 


- w n a_. sin c ,)w n sin c. + M- T (a-, sin c, 

Ci i 0 3 * X i 
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mm mm*™**?' 


i 


i 


(49c) 


+ Wg a^ cos cj w Q cos Cj]/M + ka j + Caj = 0 
j= 1,2, 3, 4 equations 

•* ** • 

Mgb. + MjWg bj + 2 w o 9 2 ) sin + c j J 

+ (-9 2 + 2Wg 0 X ) cos (« + Cj) + kbj + Cbj * 0 (49d) 

j = 1,2, 3, 4 equations 

The homogenous system of linear differential, equations can 
be expressed in state vector form:® 

X(t) = ACt + T)X(t) (50) 

where X(t) is a (20 x 1) column matrix of the coordinates and 
velocities and A(t + T) represents a (20 x 20). matrix. Typical 
periodic coefficients in A(t + T) are: 

sin(^ + c^) and cos(<|» + c^) 

which occur in Eqs. (49a), (49b) and (49d) for constant length. 

A. Application of the Lyapunov Reductibility Theorem 

The Lyapunov reductibility theorem® may be applied to reduce 
the linear system of differential equations with periodic coefficients 
to a kinematically similar system of differential equations. (A 
recent application of this theorem, in studying the stability of a 
dual-spin spacecraft with a flexible momentum wheel is given in 



Ref. 6.) The theorem states that there is a bounded nonsingular 
transformation matrix SCt) with a bounded inverse S“^Ct) such that 
a matrix D, defined by: 



D = S" 1 AS - S' 1 S 

cs: 

is a constant matrix. 

Me can then relate 



X = SCt)Y 

cs 

and obtain the system 

• 

Y = DY 

(5 


as tiie ki nemat ically equivalent of Eq. (£0). T?« bounded trans- 
formation of the form of Eq. (52) for the present system can be 

expressed: 
r -r 


9 1 

0 , 


e l 


cos ij> -sin ill 
sin ip cos ip 
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r 



In Eq. (5*0 if Z s Ce^.. then 

Y = 

Thus, the kinematically equivalent system can be expressed 
in terns of the Z matrix and the constant coefficient matrices, 

M, 6, D, K, as: 

MZ + GZ + DZ + KZ=0 (55) 

M is tite symmetric mass inertia matrix, G is a skew symmetric 
matrix involving gyroscopic terms, D is a symmetric positive 
semidefinite damping matrix and K is a symmetric stiffness matrix. 
These matrices can be expressed as follows : 
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GOO 
0 G 0 
0 0 A 
0 0 B 
0 0 -C 
0 0 -B 
D -E 0 
EDO 
-DEO 
-E -D 0 


0 0 0 D 
0 0 0 -E 
B -C -B 0 
A -B C 0 
-B A B 0 
C B A 0 
0 0 0 F 
OOOO 
0 0 0 0 
0 0 0 0 


E -D -E 

D E -D 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

F 0 0 

0 F 0 

OOF 


(58) 


( 59 ) 


where, ^ £q. (rg) j 

A = k + r Cl - Mj (sin 2 c^ - cos 2 c^J/rM] 


B = 2(H^ Wq/M) sin cos c^ 

C = Wq/M (sin 2 c ^ - cos 2 c^ 
2 

D = Wg sin c ^ 

2 

£ = % w o 003 °x 

F = ^ w 2 + k 


X 3 W 0 “ I 1 W 0 


B. Application of the Kalvin-Tait-Chetaev Theorem 
In the equation: 

MZ +DZ + GZ + KZ = 0 (60) 

we assume that M and D are synrcetric positive definite, G is 
skew symmetric, and that none of the eigenvalues of the symmetric 
matrix, K, is zero. Then the stability of the solutions of Eq. C60) 

is the same as the stability of the solutions of the truncated 

. 7 

equation 

MZ + KZ = 0. 


If M, K, and D are positive definite, the zero solution of Eq. (55) 
is asymptotically stable; if M and K are positive definite and D is 
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positive seraidefinite with (D i 0), the solution of Eq. (55) is 

7 

stable in the sense of Lyapunov but not asymptotically stable. 

The Kelvin-Tait-Chetaev Theorem can be applied to this particular 
system, where D is positive semidefinite in the absence of damping 
on the main hub, by examining the. M. and K matrices for positive defi- 
niteness. If, for a particular set of system parameters, any of 
the principal minor subdeterminants of either the M or K matrices, 
Eq. (56) cr Eq. (59), is negative, then the system will be unstable 
in the Lyapunov sense. 

C. Extension of Stability 

A recent extension of the Kelvin-Tait-Chetaev theorem to 

establish asymptotic stability of linear systems where the damping 

8 

matrix is only positive semidefinite, was presented by Muller 
and is summarized briefly, here. 

The mechanical system, 

MZ + (D + G)Z + KZ = 0 

with the fxf matrices 

M = M T > 0, D = D T > 0, G = -G T , K = K T 
is then asymptotically stable if the fxf matrix K is; (a) positive 
definite, and (b), that .the rank of an augmented matrix, S be : 

Rank of S = Rank [r*Fr • F 2 r : . . . if 1 ' 1 r ] = n (n = 2f ) 

• * • » 
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where F is a (2f x 2f) matrix and r is a (2f x f) matrix. Condition 
Ca) is a result of applying Sylvester’s theorem when the principal 
minors of matrix K are positive definite. * This extension 
of the Kelvin-Tait-Chstaev theorem which tests the system for 
asymptotic stability can be applied to this particular system. K 
can be tested for positive definiteness for a partic ular set of 
system parameter's and the augmented matrix, S, can be constructed 
as previously indicated and tested for a rank equal to n (twenty, 
for this system) . If the rank of S is n then the system will be 
asymptotically stable in the Lyapunov sense, even when there is 
no damping on the main (hub) part of the satellite. 
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CHAPTER IV 


Concluding Comments 

The rotational equations of motion for a spinning spacecraft 
system with deployable appendages and the criteria for the stability 
of such a system when the appendages are?af a fixed length, have 
been developed. The future analysis of this system will include 
the following: 

1. an attempt to identify mode shapes for constant 
extension rates, 

2. an application of the Kelvin-Tait-Chetaev theorem 
and its extension tc predict the stability of this 
system before and after deployment maneuvers for 
different values of system parameters. 

3. a computer simulation of the equations of motion 

with first order flexibility present and then 

• • 

absent, (a^ = b^ = = tu - 0) 

4. a computer simulation of the equations of motion 
with first order flexibility, in-plane bending only. 

5. an attempt to relate the El (flexural rigidity) of 
the wires tc an effective restoring linear spring 
constant, k. 
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Conclusions regarding the effect of flexibility on the system 
during extension or retraction as well as oh the stability of the 
system before or immediately after extension or retraction, 
(numerical results) will be forthcoming. 
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Fig. 2. Variables specifying position of 


wire elacent . 



